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AM  ADAPTIVE  NONCOHERENT  PROCESSOR  FOR  MJLTI -CHANNEL  RADARS 


by 

Robert  U.  Herring 
Ross  M.  Turner 
Eric  K.L.  Hung 


ABSTRACT 


An  adaptive  noncoherent  processor  (ANCP)  system,  intend¬ 
ed  as  an  ECCM  for  multi-channel  (e.g.,  frequency  agile)  radars 
operating  in  a  jamming  environment,  is  described.  The  ANCP 
system  operates  by  estimating  the  noise  power  in  each  active 
channel  of  the  radar,  adaptively  adjusting  the  channel  gain, 
and  then  noncoherently  combining  the  outputs  of  the  channels 
before  executing  a  single  target  present /not-present  threshold 
detection. 

An  example  of  the  theoretical  performance  of  an  ANCP 
system  in  the  presence  of  white-noise  jamming  with  Gaussian 
statistics  is  presented.  It  is  shown  that  a  practical  system 
performs  within  3  dB  or  better  of  an  ideal  optimum  system 
using  prior  knowlege  of  the  jamming  and  target  powers.  The 
report  concludes  with  appendices  outlining  the  mathematical 
techniques  for  computing  these  results. 


1 .  INTRODUCTION 


Multi-channel  radars,  such  as  frequency-agile  radars  or  netted 
radars,  are  used  for  their  enhanced  target-detection  and  ECCM  capabili¬ 
ties.  The  purpose  of  this  report  is  to  describe  the  theoretical  target- 
detection  and  anti-jam  properties  of  an  adaptive  noncoherent  processor 
(ANCP)  which  produces  for  each  range-doppler  cell  a  single  target 
present/non-present  decision  base  on  the  combined  outputs  of  the  multiple 
channels.  It  is  also  shown  how  to  estimate  these  properties  analytically. 

The  ANCP  operates  by  estimating  the  noise  power  present  in  each 
channel  and  then  using  this  estimated  noise  level  to  adjust  the  gain  of 
the  channel  according  to  some  predetermined  algorithm.  This  approach  is 
particularly  useful  when  the  noise  levels  differ  in  the  various  channels. 
Such  is  the  case,  for  example,  when  a  frequency-agile  radar  is  subject  to 
wlde-band  barrage  jamming  in  its  side lobes,  since  the  frequency  dependence 
of  the  detailed  sldelobe  structure  can  induce  significant  channel-to- 
channel  variation  in  the  observed  jammer  noise. 
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In  this  report  it  is  assumed  that  the  target  statistics  on  a  chan- 
nel-to-channel  basis  are  Rayleigh  or  Swerling'a  Case  2,  and  that  the  noise 
in  each  channel  is  white  with  Gaussian  statistics.  The  adaptive  gain 
algorithms  are  derived  under  these  assumptions.  It  turns  out  that  for 
practical  implementation  the  ANCP  reduces  to  a  form  of  the  ratio  detector 
described  by  Trunk  [ 1—3 j .  The  main  differences  between  the  present  work 
and  that  of  Trunk  are  in  the  derivation  of  and  justification  for  the 
choice  of  the  channel- to-channel  gain  adjustment  algorithm,  and  in  the 
techniques  used  to  estimate  the  target  detection  and  anti-jam  performance. 

Section  2  contains  a  description  of  the  operation  of  a  generalized 
ANCP  system  in  the  presence  of  white  noise  jamming.  It  begins  with  the 
derivation  of  the  optimum  adaptive  gain  function,  where  it  is  assumed  that 
prior  knowledge  of  the  average  noise  levels  is  available.  It  is  shown 
that  considerable  simplification  in  implementation  occurs  if  the  ANCP  sys¬ 
tem  is  optimized  for  large-amplitude  echoes,  in  which  case  it  becomes  a 
ratio  detector.  Section  3  describes  and  compares  the  theretically  expect¬ 
ed  performance  of  three  ANCP  systems.  The  first  is  the  optimum  ANCP  sys¬ 
tem  with  prior  knowledge  of  its  jamming  environment.  The  second  is  a 
simplified  sub-optimum  ANCP  system  with  prior  knowledge  of  its  jamming 
environment.  The  third  is  a  sub-optimum  ANCP  system  which  adapts  itself 
to  the  observed  jamming  environment.  Only  this  third  system  is  easily 
realized;  the  first  two  serve  as  benchmarks  against  which  to  measure  the 
cost  of  simplifying  the  system's  operation  and  the  cost  of  estimating  the 
noise  environment  directly.  Section  4  summarizes  these  results  and  makes 
recommendations  for  further  investigations.  This  report  also  includes  an 
Appendix  outlining  the  mathematical  procedures  followed  in  computing  the 
results  of  Section  3. 

2.  SYSTEM  DESCRIPTION 

Fig.  1  shows  an  ANCP  system  suitable  for  a  step-scan  frequency-agile 
radar  with  M  channels.  (For  clarity,  the  prcessing  for  only  a  single 
range  cell  has  been  included).  The  system  operates  as  follows:  The  sig¬ 
nal  in  each  carrier-doppler  channel  (m,k)  is  square-law  detected  to 
yield  the  sample  variate  t^(k)  defined  by  Eqn.  (1): 

M 

tM(k)  -  l  G(M,k)  t (m,k)  .  (1) 

m«l 

tf^(k)  is  then  compared  with  an  adaptively  computed  threshold,  Tjj(k),  to 
decide  whether  or  not  a  target  is  present  in  the  k'th  doppler  channel  of 
the  test  range  cell.  (The  circumflex  denotes  the  sample  value  of  a  stoch¬ 
astic  quantity.) 

2.1  Gain  Estimation 

In  Appendix  A  it  is  shown  that  for  white  noise  jamming  with  Gaussian 
statistics  and  Swerllng's  Model  2  (Rayleigh)  channel-to-channel  target 
statistics,  the  optimum  value  for  the  adaptive  gain  for  channel  (m,k), 
y0pt(m»k)»  18  given  by, 


TARGET 

/ 


Figure  1  -  Frequency- Agile  Implementation  of  an  AN CP  System. 
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Gopt(m,k) - ^ -  (2) 

[1+X(m,k)] [1+X(m,k)+S0] 

where  unit  internal  white  noise  power  spectral  density  (PSD)  has  been  * 

assumed,  X(m,k)  is  the  external  white  noise  PSD  due  to  jamming  and  is 

assumed  to  be  known  exactly,  and  SQ  is  the  particular  average  target 

signal  energy  to  quiescent  internal  noise  PSD  ratio  (SNR)  to  which  the 

system  has  been  matched. 

One  purpose  of  this  report  is  to  investigate  the  effect  of  replacing 
[1+X(m,k)]  in  Gqn.  (2)  with  estimated  values  based  on  observed  samples  of 
the  total  internal  plus  external  noise.  An  unbiased  estimate  for 
[1+X(m,k)]  is  given  by: 

1  L 

YL(m,k)  -  —  l  y£<m,k)  (3) 

L  1-1 

A 

where  the  y£(mk)*s  are  square-law  detected  noise  samples  drawn  from 
nearby  range  bins  which  are  assumed  not  to  contain  target  signals,  and  the 
subscript  L  denotes  an  L-sample  estimate.  Under  the  assumption  of 
Gaussian  statistics,  Y^(m,k)  has  a  Chi-squared  distribution  with  2L 
degrees  of  freedom.  A  suitable  system  for  implementing  Eqn.  (3)  is  shown 
in  Fig.  2. 

An  estimate  of  GopC(m,k)  is  given  by:  « 

Gopt(m»k>  "  7 - — “T2 - 

YL(m,k)  [YL(m,k)+S0] 

Equation  (4)  can  be  substituted  into  Eqn.  (1)  to  define  the  operation  of 
an  optimum  ANCP  system. 


2.2  Threshold  Estimation 

For  any  adaptive  detection  scheme  it  is  necessary  to  set  the 
detection  threshold  on  the  basis  of  the  observed  noise  and  interference 
levels.  It  will  now  be  shown  how  to  do  this  for  an  ANCP  system. 

In  general,  the  threshold  level  is  set  on  the  basis  of  providing  an 
acceptable  rate  of  false-alarm  threshold  crossings  when  there  is  known  to 
be  no  target  present.  Ideally,  this  threshold  level  can  be  estimated  as  a 
simple  function  of  the  observed  noise^levels .  The  simplest  possiDle 
estimator  for  the  adaptive  threshold  Tjj(k)  for  an  ANCP  system  would  be 
given  by: 


fM(k)  -  YZM(k) 


(5) 


L  NOISE  SAMPLES 
2  CHANNEL  k 


Figure  2  -  Noise  Monitor,  Doppler  Channel  k 


where 

M  _ 

Z(i(k)  -  l  G(m,k)  YL(m,k)  (6) 

m-1 

and  y  is  a  constant  to  be  determined  on  the  basis  of  the  statistical 
distribution  of  Zfl(k).  Zfl(k)  is  an  estimate  of  the  adaptively  combined 
noise  powers  of  the  threshold  detector. 

A  fundamental  difficulty  in  using  the  estimated  optimum  gain 
function  in  practice  is  revealed  if  Eqn.  (6)  is  further  analyzed  by 
considering  the  contribution,  z(m,k),  of  a  single  channel  to  the  sum  in 
Eqn.  (6),  where 


If  Che  optimum  gain  function  of  Eqn.  (4)  is  substituted  into  Eqn.  (7), 
then, 

z0pt(®.k)  ”  *-~+Sq -  (8) 

YL(m,k)  +  S0 

and  it  is  clear  that  Z^opt^K)  as  defined  by  Eqns.  (7)  and  (6)  is  the 
sum  of  M  non-linearly  transformed  Chi-squared  distributed  variates,  each 
with  2L  degrees  of  freedom  and  potentially  different  means. 

It  seems  apparent  that  the  statistical  distribution  of  Z ^  0pC(k) 
is  not  easily  calculable,  nor  is  it  likely  that  the  desirable  simple 
linear  relationship  between  the  threshold  setting  and  the  noise-only 
output  variate  as  defined  by  Eqn.  (5)  is  valid.  However,  this  challenging 
problem  is  obviated  by  optimizing  for  very  large  SNRs,  which  results  in  an 
easily  implemented  system  with  sub-optimum  performance. 

In  the  limit  SQ  *  “,  Eqn.  (8)  reduces  to 

ZoaCm.k)  ■  1  (9) 

so  that,  according  to  Eqn.  (5),  the  threshold  T(j,»(k)  is  given  by 

TM,„.(k)  -  yM  (10) 

which  is  a  constant  Independent  of  the  observed  external  noise  level. 

This  result  means  that  in  principle  an  ANCP  system  using  inverse  noise 
weights  gives  constant  false-alarm  rate  (CFAR)  performance  using  a  fixed 
threshold,  because  of  the  normalization  of  the  M  components  of  the  test 
variate  t^(k).  In  practice,  this  result  means  that  the  threshold  can  be 
set  empirically  and  can  be  expected  to  be  relatively  stable,  provided  the 
assumption  of  Gaussian  noise  statistics  is  not  badly  violated.  This 
remains  as  a  point  for  experimental  investigation. 

Again  in  the  limit  S0+®,  the  gain  function  Goo(m,k)  is  found  from 
Eqn.  (4)  to  be: 


G«(m,k) 


YL(m,k) 


(ID 


which  is  simple  to  estimate  and  implement.  For  these  compelling  reasons, 


Che  remainder  of  this  report  shall  be  concerned  principally  wich  Che  per¬ 
formance  of  Che  sub-optimum  inverse-noise  weighted  implementation  of  the 
ANCP  system,  although  results  for  the  optimum  system  have  been  included  as 
the  ultimate  standard  of  performance. 

3.  SYSTEM  PERFORMANCE 

In  this  section  the  theoretically  expected  detection  and  anti-jam 
performance  of  an  ANCP  system  is  investigated.  The  number  of  carrier  fre¬ 
quencies,  M,  has  arbitrarily  been  chosen  to  be  10. 

Section  3.1  contains  discussions  of  the  theoretical  performance  of 
ANCP  systems  when  the  external  noise  powers  are  known  exactly.  Section 
3.1.1  discusses  the  performance  of  a  system  using  the  optimum  gain  func¬ 
tion.  These  discussions  provide  a  basis  for  quantitative  estimates  of  how 
much  the  system  performance  is  degraded  by  the  use  of  the  sub-optimum 
inverse  gain  function. 

Section  3.2  contains  results  of  the  theoretical  analysis  of  the  per¬ 
formance  of  10-channel  ANCP  systems  using  Inverse  gains  estimated  using 
L*8,  16  and  32  samples  of  the  external  noise  in  each  of  the  channels. 
Comparison  of  these  results  with  those  of  Sections  3.1.1  and  3.1.2  allows 
estimation  of  the  losses  due  to  using  sample  observations  of  the  ambient 
noise  instead  of  prior  knowledge. 

In  both  Sections  3.1  and  3.2,  the  numerical  results  were  computed 
using  the  characteristic  function  method  of  Bird  [3].  In  Section  3.1  the 
required  characteristic  functions  are  straightforward  modifications  of  the 
well-known  characteristic  function  of  the  Chi-squared  distribution.  In 
Section  3.2  matters  were  considerably  more  complicated.  The  mathematical 
details  are  contained  in  Appendix  B. 

In  the  following  analyses  these  basic  assumptions  were  made:  The 
target  echo  amplitudes  were  assumed  to  be  uncorrelated  on  a  channel-to- 
channel  basis,  and  the  target  echo  amplitudes  were  assumed  to  have 
Swerling's  Case  2  or  Rayleigh  statistics.  It  was  assumed  that  the  powers 
radiated  at  each  carrier  frequency  were  identical,  so  that  the  received 
echo  decorrelation  on  a  channel-to-channel  basis  was  due  solely  to  the 
target  cross-section  characteristics.  (The  effect  of  a  systematic 
channel-to  channel  variation  in  target  echo  strength  can  be  incorporated 
into  the  analysis,  as  pointed  out  in  Appendix  B.3.1.  Such  variations 
arise,  for  example,  due  to  scanning  modulation  in  rotating  antenna 
radars . ) 

The  system  chosen  for  analysis  had  10  channels  (M**10)  with  a  design 
false-alarm  rate  (FAR)  of  10”6  and  a  probability  of  target  detection 
(Pq)  of  90  percent.  The  number  of  jammed  channels,  J,  ranged  from  0  to 
10,  and  all  the  jammed  channels  were  assumed  to  have  the  same  enhanced 
noise  power  levels.  The  parameter  chosen  as  a  measure  of  system  perform¬ 
ance  was  S,  the  average  received  echo  energy  in  a  single  channel.  S  was 
expressed  in  decibels  above  the  unit  internal  white  noise  PSD  of  an  un¬ 
jammed  channel. 
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3.1  Performance  of  ANCP  Systems  In  a  Known  Noise  Environment 

3.1.1  Optimum  Gain  Function 

Results  for  an  ANCP  system  with  perfect  prior  knowledge  of  its  noise 
environment  are  shown  in  Fig.  3.  The  number  of  jammed  channels,  J,  has 
been  plotted  along  the  horizontal  axis  and  S,  the  single-channel  average 
echo  energy  acquired  to  produce  a  Pq  of  90  percent  where  the  threshold 
has  been  set  to  maintain  a  FAR  of  10“  ,  has  been  plotted  along  the  verti¬ 
cal  axis.  The  system  performance  has  been  plotted  as  curves  of  S  vs  J  for 
levels  of  internal  plus  external  noise  in  the  jammed  channels  ranging  from 
3  dB  to  30  dB  above  the  quiescent  internal  noise  level  of  0  dB.  At  each 
point  on  these  curves,  the  value  of  S0  in  the  optimum  gain  function 
Gopt  of  Eqn.  (2)  was  set  equal  to  S.  The  values  of  S  required  by  sys¬ 
tems  composed  of  M**l  to  M”10  channels  to  maintain  this  same  performance  in 
the  absence  of  jamming  have  also  been  indicated  in  this  figure. 


L  *r*« 


V 

* 


& 


,Vv 


r-ir. 


In  Fig.  3  (and  in  the  other  performance  curves  to  be  discussed 
below),  there  are  certain  fixed  points  of  Interest.  Specifically,  these 
points  correspond  to  the  case  of  no  jamming  (J=0)  and  all  channels  jammed 
equally  (J**10).  In  each  of  these  cases  the  gains  in  all  10  channels  be¬ 
come  identical,  regardless  of  the  gain-setting  algorithm  used.  (In  the 
case  of  all  channels  jammed  identically,  S  mist  be  increased  by  the  same 
number  of  dB  as  the  noise  level,  in  order  to  maintain  the  required  P^.) 

The  effect  of  the  choice  of  adaptive  gain  function  is  to  control  the  shape 
of  the  curves  joining  these  fixed  points. 


The  important  observations  from  Fig.  3  are  that  in  the  presence  of 
jamming  the  optimum  ANCP  system  performs  at  least  as  well  as  one  having 
10-J  channels,  and  that  the  J  jammed  channels  are  more  or  less  Ignored  by 
the  system,  according  to  how  heavily  they  are  jammed.  In  the  case  of 
"light”  jamming  (10  dB  or  less),  the  jammed  channels  continue  to  make 
a  useful  contribution  to  the  summed  output.  In  the  case  of  "heavy"  jam¬ 
ming  (20  dB  or  more),  the  jammed  channels  are  effectively  "turned  off"  and 
ignored. 


3.1.2  Inverse  Noise-Power  Gain  Function 


Results  for  an  ANCP  system  which  has  perfect  prior  knowledge  of  its 
noise  environment  but  which  has  been  optimized  for  a  fixed  value  of  SQ  51 
«  have  been  plotted  in  Fig.  4.  It  was  shown  in  Section  2.2  that  such  a 
system  has  its  channel  gains  set  inversely  proportional  to  the  white-noise 
PSD  in  each  channel.  The  required  system  performance  of  FAR 


1U“6  and 


90  percent  is  the  same  as  that  of  Fig.  3. 


Comparison  of  Figs.  3  and  4  shows  the  following.  First,  the  values 
for  S  at  the  fixed  points  (J*0  and  J“10)  are  the  same  for  both  figures,  as 
discussed  in  Section  3.1.1.  Second,  for  light  jamming  (1U  dB  or  less), 
the  inverse  gain  function  gives  performance  within  0.4  dB  of  that  of  the 
optimum  gain  function. 


FAR  «  10< 

Pq  *  90% 

CASE  2 

10  CHANNELS 

INVERSE  weighting 

EXACT  KNOWLEDGE  OF  THE  NOISE  (L=»l 


NUMBER  OF  NO'SV  CHANNELS.  J 


Figure  4  -  Required  Average  Single-Channel  Echo  Energy  for  a  10-Channel 
AN CP  Systea  Using  a  Priori  Knowledge  of  the  Noise  Environnent 
and  fhw  Inverse  Noise  Gain  Function. 


For  heavy  jamming  (20dB  or  greater)  the  performance  of  the  inverse 
gain  system  is  inferior  by  amounts  ranging  from  0.6  dB  for  J*5  to  2  dB  for 
J~9.  The  reason  for  this  is  that  under  heavy  jamming  the  outputs  from  the 
jammed  channels  are  not  suppressed  or  “turned  off",  but  instead  each  chan¬ 
nel  contributes  an  average  of  one  unit  of  noise  to  the  accumulator,  there¬ 
by  degrading  the  SNR  at  the  threshold  detector.  However,  It  is  the  pres¬ 
ence  of  this  constant  noise  power  level  at  the  threshold  detector  that 
leads  to  the  simple  fixed-threshold  CFAR  performance  of  this  system.  For 
the  10-channel  system  considered  here,  this  CFAR  performance  costs  2  dB  or 
less  in  terms  of  additional  power  required. 

3.2  Performance  of  ANCP  Systems  which  Estimate  their  Noise  Environment 

Results  for  an  ANCP  system  which  uses  inverse  noise  gains  and  esti¬ 
mates  its  noise  environment  according  to  Eqn.  (3)  and  an  implementation  of 
the  scheme  diagrammed  in  Fig.  2  are  given  in  Fig.  5.  The  value  of  L,  the 
number  of  noise  observations  per  channel,  was  chosen  to  be  16.  All  other 
parameters  were  as  before. 

Comparison  of  Figs.  4  and  3  shows  that  the  shapes  of  the  character¬ 
istic  curves  are  essentially  identical,  but  that  those  of  Fig.  3  are  dis¬ 
placed  upwards  by  a  factor  of  about  1.0  dB.  This  factor  of  1.0  dB  is  the 
loss  incurred  by  not  having  exact  knowlege  of  the  white  noise  environment. 

The  cases  for  L-8  and  L-32  were  also  investigated.  For  the  case 
L-32,  the  estimation  loss  was  about  0.4  dB;  for  the  case  L»8  it  was  about 
3  dB.  The  implication  is  that  for  a  10-channel  system,  16  noise  samples 
per  channel  provide  adequate  information  about  the  noise  environment, 
while  at  the  same  time  allowing  rapid  adaptation  in  a  non-stationary  noise 
environment. 


4.0  SUMMARY  AMD  RECOMMENDATIONS 

An  adaptive  non-coherent  processor  (ANCP)  system  for  use  in  multi¬ 
channel  radar  ECCM  applications  has  been  described  and  its  performance  has 
been  analyzed  theoretically.  Of  the  three  forms  of  Implementation  ana¬ 
lyzed,  the  first  two  took  advantage  of  prior  knowledge  of  the  noise  and 
jamming  environment  and  thus  are  of  academic  interest  only.  However,  they 
provided  benchmark  information  in  assessing  the  performance  of  the  third 
implementation,  which  uses  observations  of  the  noise  environment  to  adapt 
its  operating  characteristics.  It  was  shown  that  for  the  particular  sys¬ 
tem  investigated,  having  10  Independent  frequency  channels,  a  false-  alarm 
rate  of  10"  and  a  target  detection  probability  of  90  percent,  the  in¬ 
crease  in  radiated  power  needed  to  compensate  for  both  a  non-optimum 
operating  characteristic  or  adaptive  gain  function  and  a  lack  of  prior 
knowledge  of  the  jamming  and  noise  environment  was  a  modest  3  dB. 

It  is  recommended  that  further  work  relevant  to  the  ANCP  systm  be 
done  in  three  general  areas.  The  first  is  to  compare  the  performance  of 
the  ANCP  system  with  the  performance  of  systems  which  perform  threshold 
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Figure  5  -  Required  Average  Single-Channel  Echo  Energy  for  a  10-Channel 
AMCP  Systea  Using  16  Saaples  of  Che  Noise  EnvlronaenC  in  Each 
Channel  to  Estimate  the  Inverse  Noise  Gains. 


detection  on  their  individual  channels  and  then  process  the  detector  out¬ 
puts  in  some  manner.  One  such  candidate  system,  which  has  already  been 
investigated  [4],  is  the  M-channel  system  having  a  conventional  adaptive- 
threshold  CFAR  detector  in  each  of  the  M  channels.  The  adaptive  thresh¬ 
olds  are  set  to  maintain  a  desired  output  false-alarm  rate  and  a  target 
detection  is  declared  when  the  threshold  levels  are  exceeded  in  N  or  more 
channels.  N  must  be  selected  to  optimize  system  performance  in  some 
manner. 

The  second  general  area  of  investigation  is  to  simulate  the  perform¬ 
ance  of  an  ANCP  systm  via  software  modelling.  In  this  manner  any  undue 
sensitivity  of  system  performance  to  the  presence  of  forms  of  jamming 
other  than  white  noise,  such  as,  for  example,  blinking  tone  jammers  or 
frequency-sweeping  jammers,  would  be  revealed. 

The  third  general  area  of  investigation,  subject  to  favourable  re¬ 
sults  being  obtained  for  the  first  two,  is  to  perform  measurements  using 
an  actual  radar  system  such  as  the  CKC  phased-array  radar.  Successful 
results  for  these  experiments  would  indicate  that  development  of  an  add-on 
processor  for  existing  frequency-agile  radar  systems  would  be  warranted. 
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APPENDIX  A 


i 

DERIVATION  OF  TUB  OPTIMUM  GAIN  FUNCTION  j 


For  Swerling's  Case  2,  Che  echoes  and  the  noises  are  mutually 
uncorrelaced  within  a  channel,  and  the  channels  are  also  uncorrelated. 
Thus,  for  a  two-channel  ANCP,  the  covariance  matrix  for  the  k'th  doppler 
channel  is 


r 


[Rl  ] 


1+X( 1 ,k)+S0 


1+X(2,k)+S0 


(A.l ) 


when  a  signal  with  energy  S0  is  present.  If  Che  signal  is  not  present, 
then  the  covariance  matrix  is 


[Bq] 


” 1+X( 1 ,k)  0 

0  1+X(2,k) 


(A. 2) 


Let  [i(i]  ■  [Ro]_1~lRi ]_1  be  the  difference  between  the  Inverses 
of  [Ry]  and  [R^J.  Then,  since  both  the  echoes  and  noises  have  been 
assumed  to  have  Gaussian  amplitude  distributions,  the  test  statistic  for 
making  a  Neyman-Pearson  decision  regarding  signal  presence  or  absence  in 
the  k'th  doppler  channel  is 


Tk  *  r£>  1*1  ^  1  J 

where  rj  and  r£  are  the  outputs  of  the  ANCP  doppler  channels  (l,k)  and 
(2,k).  Eqn.  (A. 3)  can  be  rewritten  as 

2  s  2 

Tk  *  l  - ^ -  | rm|  2  CA.4) 

m+1  [ 1+X(m,k) ] [ 1+X(m,k)+SQ  ' 

Extension  to  more  than  two  doppler  channels  is  trivial. 

A  simple  scaling  change,  replacing  SQ  by  (1+SQ)  in  the  numerator 
of  the  fraction  in  Eqn.  (A.4),  gives  Eqn.  (2)  for  the  optimum  gain 
(*opt(|D*lc)*  only  effect  of  the  scale  change  is  to  give  G0pt(m,k) 

unit  value  when  there  is  no  external  noise.  * 
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APPENDIX  B 

MATHEMATICAL  DETAILS 


B.l  INTRODUCTION 

The  numerical  results  of  Section  3  were  computed  using  the  method  of 
Bird  [5].  This  is  a  technique  by  which  the  cumulative  probability  distri¬ 
bution  of  the  test  statistic  at  the  accumulator  output  can  be  obtained 
directly  using  a  Fourier  series  and  the  characteristic  function  of  the 
probability  density  distribution. 

The  purpose  of  this  Appendix  is  to  show  the  derivations  of  the 
required  characteristic  functions  for  the  cases  of  prior  Icnowledge  and  of 
direct  observation  of  the  noise  environment,  and  to  show  how  the  two  error 
parameters  in  Bird's  technique  are  determined. 

B.2  DEFINITIONS  OF  FUNDAMENTAL  RESULTS 


In  order  to  proceed  it  is  necessary  to  introduce  some  definitions 
and  fundamental  results.  To  begin,  if  p(x)  is  a  probability  density 
distribution,  then  its  characteristic  function,  C(a>),  is  defined  as  its 
Fourier  transform: 


C(oj)  -  /  p(x)  e~^ UXdx  (B.l) 

o 

The  cumulative  probability  distribution  P(T),  i.e.  the  probability 
PR.{x<T},  is  defined  as: 

P(T)  -  /T  p(x)  dx  (B.2) 

o 

so  that  if  p(x)  corresponds  to  the  distribution  of  the  test  variate  at  a 
threshold  detector  and  T  is  the  threshold,  then  P(T)  is  the  probability  of 
correct  dismissal  for  the  case  of  target  not  present  or  the  probability  of 
missed  detection  for  the  case  of  target  present. 

Suppose  now  that  a  variable  gain  is  introduced  into  the  system,  so 
that  p(x)  is  replaced  by  a  conditional  probability  density  distribution 
p(yjG),  where 

y  ■  Gx  .  (B.3) 

Further,  let  it  be  assumed  that  G  is  a  stochastic  variable  with  probabil¬ 
ity  density  distribution  p(G),  and  that  x  and  G  are  statistically  indepen¬ 
dent.  It  is  desired  to  compute  the  characteristic  function  of  the  proba¬ 
bility  density  distribution  underlying  P(T),  where  now 
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P(T) 


•  T 

/  p(G)  /  p(y|G)  dy  dG 
o  o  1 


(B.4a) 


-  /  p(G)  P(T|G)  dG  (B.4b) 

The  required  probability  density  distribution  is  found  by  differentiating 
Gqn.  (B.4b)  with  respect  to  T: 


p(x)  -  /  p(G)  p(x|G)  dG  (B.5) 

so  that  the  required  characteristic  function  can  be  found  by  applying 
Eqn.  (B.l): 

ZT(oi)  -  /  p(G)  C(ujG)  dG.  (B.6) 

The  overbar  here  denotes  averaging  over  the  gain  distribution  and 

C(u>  |G)  -  £  p(x|G)  e~J ^dx  (B.7) 

is  the  characteristic  function  corresponding  to  P(xjG).  This  result  is 
used  in  Sections  B.2.1  and  B.3.1. 

B.2  Perfect  Prior  Knowledge  of  the  Noise  Environment 
B.2.1  Derivation  of  the  Characteristic  Function 


For  the  single-hit  or  single-channel  case  with  Swerling  Case  2  stat¬ 
istics,  che  characteristic  function  of  the  probability  density  function  at 
the  output  of  the  square-law  detector  is  given  by: 

C^( uj)  -  - - -  (B.8) 

l+ju>ll+S] 

where  S  is  the  average  target  SNR  and  the  average  noise  power  has  been 
normalized  to  unity  ([b],  Eqn.  (111.8)). 

If  additional  white  noise,  X(m,k),  Is  added  to  the  input  signal, 
then  the  unit  noise  Is  replaced  by  1+X(m,k)  and  the  characteristic  func¬ 
tion  for  che  single  channel  case  becomes: 


_ 1 _ 

l+ju)[l+X(m,k)+S] 


(B.9) 


C^oi) 
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If  a  deterministic  gain,  G(m,k),  is  then  incorporated  into  the  chan¬ 
nel,  the  characteristic  function  is  further  modified  to  become: 

Ci(iojG)  -  - - -  (B.  10) 

1  1+j  uG(m,k) [ 1+X(m,k)+S] 

However,  since  the  gain  is  deterministic,  its  probability  density  distri¬ 
bution  is  a  Dirac  delta  function  and  from  Eqn.  (B.6), 

T?!  (to)  -  Cx  (uj  jG)  (B.ll) 

Since  the  channels  are  assumed  to  be  statistically  independent,  the 
well-known  convolution  theorem  of  statistics  can  be  invoked  [5],  so  that 
the  characteristic  function  for  the  test  statistic  composed  of  the  sum  of 
the  M  independent  outputs  of  an  M-channel  system  is  given  by: 

M  1 

Cm,k(<*>)  -  n  - 1 -  (B.  12) 

m“1  1+j wG(m,k) [ 1+X(m,k)+S] 

B.2.2  Bounding  the  Errors 

It  is  clear  from  comparison  of  Eqn.  (B.12)  with  Eqn.  (24)  of  14] 
chat  the  expansion  interval  R  for  Che  Fourier  series  used  in  Bird's  tech¬ 
nique  is,  by  analogy  with  Eqn.  (31)  of  [5 J , 

M 

R  *  ln(M/eR)  l  G(m,k) [X(m,k)+S]  (B.13) 

m-1 

where  eR  is  chosen  as  the  largest  error  acceptable  due  to  replacing  the 
normally  Infinite  upper  limit  of  integration  in  Eqn.  (B.l)  with  the  large 
but  finite  value  R.  Appropriate  selection  of  a  value  for  R  controls  the 
first  source  of  error  in  Bird's  technique.  Values  of  eR  of  FAR/ 100  and 
10~6  were  chosen  when  setting  the  threshold  and  computing  Pp,  respect¬ 
ively. 


The  second  source  of  error  arises  due  to  truncating  the  infinite 
series  which  occurs  in  the  discrete  Fourier  transform.  Rather  than  fol¬ 
lowing  the  complicated  procedure  of  [5],  the  simple  expedient  of  stopping 
when  the  magnitude  of  the  terms  being  accumulated  in  the  sum  Eqn.  (12)  of 
[3]  become  smaller  than  eR/lO  was  adopted.  Experience  has  shown  that 
this  stopping  rule  produces  reliable  results  for  Swerling's  Case  2,  both 
when  setting  the  threshold  level  T  and  when  estimating  Pp  as  a  function 
of  T  and  S. 

B.3  Direct  Observations  of  the Noise  Environment 


B.3.1  Derivation  of  the  Characteristic  Function 


The  required  characteristic  function  can  be  derived  directly  from 


A 
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Eqns.  (B.6)  and  (B.10)  and  then  applying  the  convolution  theorem  to  derive 
a  final  result  analogous  to  Eqn.  (B.12). 

To  begin,  a  simplified  notation  will  be  introduced.  In  particular, 
indices  and  subscripts  will  be  suppressed.  Thus  Eqn.  (B.10)  is  rewritten 


■C(u)jG) 


(B.14) 


1+j  uX»[X+S  J 


where  X  is  the  average  total  (internal  plus  external)  noise  in  the  test 
range  cell,  and  the  gain  G  is  a  stochastic  variable.  The  desired  averaged 

characteristic  function,  "C'C  to) ,  is  given  by 


C(u>)  -  J 


l+juiG[X+s] 


(B. 15) 


where  p(G)  is  the  probability  density  distribution  of  the  stochastic  gain 
function. 

For  the  case  of  inverse  noise  gains,  the  gain  function  is  given  by 
Eqn.  (11): 


where 


G  -  1/Y 


*  -  l  y 

L  1-1 


(B. 16) 


(B. 17 ) 


is  the  sum  of  L  Chi-squared  variates,  each  with  two  degrees  of  freedom. 
Therefore,  Y  has  a  Chi-squared  probability  density  distribution  with  2L 
degrees  of  freedom. 

Equation  (B.15)  can  be  rewritten  in  this  case  as 


C(w)  -  / 


l+ju)[X+S]/Y 


(B. 18) 


where 


(L /T)L  YL_1  e-LY/7 


V- 


«  V 


mmMtm 


is  the  density  function  and  Y  is  the  expected  value  of  Y.  If  it  is  as¬ 
sumed  that  the  noise  in  the  test  cell  is  statistically  Identical  to  the 
noise  in  the  noise-observation  cells,  then 

Y  «  Y  .  (B.20 

All  that  remains  is  to  substitute  equations  (B.19)  and  (B.20)  into 
Eqn.  (B.18)  and  to  solve  by  repeated  application  of  integration  by  parts. 
Doing  so  yields 

_  i  f  L-1  t 

C(o»)  -  — - . —  (  l  (L-l-i) !  C-jLM'r 

(L-l)!  £»o 


+  (-jU>’)L  eJ*“'  EiCjLto')}  (B.21 ) 

where 

u'  -  u(l+S/Y),  (B.22 ) 

and  in  general 

E  .  (z)  -  r  dt  (B.23) 

n-r  i  ^  £  n 

is  the  exponential  integral  of  order  n  ([7],  Eqn.  (5.1.4)). 

Equation  (B.21)  can  be  further  simplified  by  repeated  application  of 
Eqn.  (5.1.14)  of  [5],  i.e., 

En+1(Z)  -  -  [e-z  -  z  En(Z)]  (B.24) 

n 

(n  m  1,2,3,...) 

to  get 

C(to)  -  L  ejL“’  EL+1(jLu)')  (B.25) 

where  again  <*>'  is  given  by  Eqn.  (B.22). 

It  is  noted  that  the  effects  of  systematic  channel-to-channel  vari¬ 
ations  in  the  received  signal  can  be  incorporated  by  appropriate  selection 
of  the  value  for  S  in  Eqn.  (B.22)  for  each  channel  component  of  the  char¬ 
acteristic  function. 
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B.3.2  Computacion  of  the  Characteristic  Function 

Unfortunately,  the  exponential  integral  with  complex  or  purely 
Imaginary  argument  and  order  greater  than  one  seems  not  to  have  been  a 
popular  topic  of  research,  and  an  extensive  literature  seach  did  not 
locate  any  numerical  algorithms  for  its  computation.  Hence  it  became 
necessary  to  develop  some  techniques  here. 

It  was  found  experimentally  that  the  problem  could  be  broken  into 
two  parts,  according  to  whether  w'L  was  greater  than  or  less  than  18.  For 
oj'L  <18,  it  was  found  that,  by  using  double-precision  (72-bit  or  18  deci¬ 
mal  place)  floating  point  arithmetic  on  a  Honeywell  CP-6  system,  C(  o>) 
could  be  computed  to  12-digit  accuracy  directly  from  Eqn.  (B.21). 

E^(jLui')  was  computed  by  noting  that 

ejLa>'  Ex  (jLw' )  -  g(Laj')  -  jf(Lu)')  (B.26) 

where  the  auxiliary  functions  f(Z)  and  g(Z)  are  defined  as: 

00 

f (Z)  -  /  dt  (B.27) 

o  t+Z 

and 

g(Z)  -  /  ^2L_£.  dt  U. 28) 

o  t+z 

from  Eqns.  (5.2.12)  and  (5.2.13)  of  [7J.  f(Z)  and  y(Z)  were  computed 
using  a  Chebyshev  polynomial  expansion  [8J. 

For  Lid'  ^  18,  Eqn.  (B.25)  was  used.  From  [7 J ,  Eqn.  (5.1.51),  it  is 
noted  the  En(Z)  can  be  expanded  in  a  divergent  asymptotic  series: 

En(z)  SZI  (1  -  a  +  aCfl±D.  -  ttU+Uin+a)+  ...}  (B.29) 

z  z  z2  z3 

Substitution  of  Eqn.  (B.29)  into  (B.25)  and  some  algebraic  manipulation 
yields 


{  i  +  iltlL  + 

jo)'  (-j(*l*)  (-jO)' )2 

(l+L'Sd-t-ZL'Sd+aL-1)  +  ...}  (B.30) 

(-ju)')3 


as  an  asymptotic  expansion  for  C(u>). 
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It  turns  out  that  Eqn.  (B.30)  is  compuationally  useful,  in  spite  of 
its  divergence.  A  form  of  Stieltjes  summation,  known  as  the  eta-algorithm 
[9],  was  applied  to  Eqn.  (B.30)  with  a  large  measure  of  success. 
(Stieltjes  summation  is  a  technique  for  forcing  convergence  on  a  class  of 
asymptotic  expansions  which  are  known  to  correspond  to  functions  which 
have  finite  values).  It  was  found  necessary  to  apply  the  eta-algorithm 
twice  for  the  case  of  values  of  Lw*  in  the  lower  range  before  convergence 
to  12-digit  accuracy  was  achieved.  Caution  was  advised  in  [9]  regarding 
this  procedure;  however,  here  it  appeared  to  work  successfully. 

Finally,  as  a  reassuring  note  of  consistency  it  is  observed  that  in 
the  limit  L+®,  corresponding  to  perfect  knowledge  of  the  noise  environ¬ 
ment,  Eqn.  (B.30)  reduces  to: 


-r~r  I 

(B.31a) 

n-0 

(jmT1 

l+(ju>Tl 

(B.31 b) 

1 

(B.31c) 

l+jo)(Y+S)/T 

which  is  the  same  result  as  given  by  Eqn.  (B.10)  for  the  case  of  inverse 
noise  gains  and  prior  knowledge  of  the  noise  environment. 

B.3.3  Bounding  the  Errors 

In  finding  the  value  for  the  threshold  setting,  T,  needed  to  give 
the  specified  FAR,  the  value  of  R  was  initially  set  to  5T,  where  T  was  a 
trial  estimate,  and  a  trial  value  of  the  FAR  was  determined.  T  and  R  were 
then  adjusted  iteratively  until  the  specified  FAR  was  achieved.  Finally, 
stability  of  the  result  was  confirmed  by  increasing  R  by  an  arbitrary 
factor  of  2.3  and  verifying  that  this  did  not  affect  the  estimated  value 
of  the  FAR  by  a  relative  change  of  more  than  1  percent. 

In  computing  Pq,  the  initial  value  for  R  was  chosen  to  be  the 
greater  of  either  5T  or  that  value  given  by  Eqn.  (B.13)  with  ejj  *  10-5. 

R  was  then  successively  doubled  until  the  computed  value  of  Pq  remained 
stable  to  1  part  in  10  . 

In  both  situations,  setting  T  and  determining  Pq,  the  same  stop¬ 
ping  rule  for  terminating  the  summation  was  used;  i.e.  that  the  absolute 
magnitude  of  the  last  term  in  the  truncated  sum  should  be  less  than  0.1 
times  the  required  tolerance  in  the  accuracy  of  the  FAR  and  Pq  esti¬ 
mates.  Thus  these  magnitudes  were  10-9  and  10”6,  respectively,  for  the  { 

two  situations. 
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An  adaptive  noncoherent  processor  (ANCP)  system,  intended  as  an 
ECCM  for  multi-channel  (e.g.,  frequency  agile)  radars  operating  in  a 
jamming  environment,  is  described.  The  ANCP  system  operates  by  esti¬ 
mating  the  noise  power  in  each  active  channel  of  the  radar,  adaptively 
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puts  of  the  channels  before  executing  a  single  target  present/not- 
present  threshold  detection. 

An  example  of  the  theoretical  performance  of  an  ANCP  system  in 
the  presence  of  white-noise  jamming  with  Gaussian  statistics  is 
presented.  It  is  shown  that  a  practical  system  performs  within  3  dB 
or  better  of  an  ideal  optimum  system  using  prior  knowlege  of  the 
jamming  and  target  powers.  The  report  concludes  with  appendices 
outlining  the  mathematical  techniques  for  computing  these  results. 


